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Abstract. We study the magnitude- redshift relation for the Type la super- 
! novae data and the angular size-redshift relation for the updated compact 

radio sources data (from Gurvits et al) by considering four variable A-models: 
A ~ S~ 2 , A ~ H 2 , A ~ p and A ~ t~ 2 . 
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It is found that all the variable A-models, as well as the constant A- 
Friedmann model, fit the supernovae data equally well with % 2 /dof ~ 1 and 
require non-zero, positive values of A and an accelerating expansion of the 
universe. The estimates of the density parameter for the variable A-models 
are found higher than those for the constant A-Friedmann model. 

From the compact radio sources data, it is found, by assuming the no- 
evolution hypothesis, that the Gurvits et al' model (Friedmann model with 
A = 0) is not the best-fitting model for the constant A case. The best-fitting 
Friedmann model (with constant A) is found to be a low density, vacuum- 
dominated accelerating universe. The fits of this data set to the (variable, 
^ ■ as well as, constant A-) models are found very good with x 2 /dof ~ 0.5 and 

require non-zero, positive values of A with either sign of the deceleration 
parameter. However, for realistic values of the matter density parameter, 
the only interesting solutions are (a) estimated from the supernovae data: 
the best-fit solutions for the flat models (including the constant A case); (b) 
estimated from the radio sources data: the global best-fit solutions for the 
models A ~ H 2 and A ~ p, the best-fit solution for the flat model with A = 
constant and the Gurvits et al' model. 

It is noted that, as in the case of recent CMB analyses, the data sets seem 
to favour a spherical universe {k > 0). 

PACS numbers: 04.20. Jb, 98.80.Es 
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1. Introduction 

The cosmic distance measures, for example, the luminosity distance and the 
angular size distance, depend sensitively on the spatial curvature and the 
expansion dynamics of the models and consequently on the parameters of the 
models. For this reason, the magnitude-redshift (m-z) relation for a distant 
standard candle and the angular size-redshift (Q-z) relation for a distant 
standard measuring rod have been proposed as potential tests for cosmological 
models and play crucial role in determining cosmological parameters. 

The observations on magnitude and redshift of Type la supernovae made, 
independently, by Perlmutter et al (1999) and Riess et al (1998) appear to 
suggest that our universe may be accelerating with a large fraction of the 
cosmological density in the form of the cosmological A-term. Their findings 
arise from the study of more than 50 Type la supernovae with redshifts in 
the range 0.16 < z < 0.83 and suggest Friedmann models with negative 
pressure- matter such as a cosmological constant, domain walls or cosmic 
strings (Vilenkin 1985, Garnavich et al 1998). 

On the other hand, models with a dynamic cosmological term A(t) has 
been considered in numerous papers to explain the observed small value of 
A, which is about 120 orders of magnitude below the value for the vacuum 
energy density predicted by quantum field theory (Wienberg 1989, Carroll 
et al 1992, Sahni and Starobinsky 2000). It has been argued that, due to 
the coupling of the dynamic degree of freedom with the matter fields of 
the universe, A relaxes to its present small value through the expansion 
of the universe and the creation of photons. This approach is essentially 
phenomenological in nature but explains, in a natural way, the present small 
value of A which might be large in the early universe. From this point of 
view, the cosmological constant is small because the universe is old. 

As the dynamics of the variable A-models depends sensitively on the 
chosen dynamic law for the variation of A and, in general, becomes altogether 
different from the dynamics of the corresponding constant A-models, there is 
no reason to believe that the observations of distant objects would also agree 
with the variable A-models, given that they agree with the corresponding 
constant A-ones, especially for the same estimates of the parameters. In this 
view, it would be worth while to test the consistency of these observations 
with the variable A-models and find the estimates of different cosmological 
parameters required by these models. 

In this paper, we consider the data on the magnitude and redshift of 
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Type la supernovae from Perlmutter et al and the data on the angular size 
and redshift of compact radio sources updated and extended recently by 
Gurvits et al (1999) and study the m-z and the Q-z relations in some decaying 
A-models. For this purpose, we consider the following four cases of the 
phenomenological decay of A: 

case (1): A ~ S~ 2 , 

case (2): A ~ H 2 , 

case (3): A ~ p, 

case (4): A ~ t~ 2 , 

where S, H, p and t are, respectively, the scale factor of the R-W metric, 
the Hubble parameter, the energy density and the cosmic time. These are 
among the main dynamical laws one finds in the literature proposed for the 
decay of A. Case (1), which has been the subject of most attention (Ozer 
and Taha 1987; Chen and Wu 1990; Abdel-Rahman 1992; Abdussattar and 
Vishwakarma 1996, 1997; to mention a few), was originally proposed by 
Chen and Wu through dimensional arguments made in the spirit of quantum 
cosmology. It is interesting that this ansatz also comes out as a natural 
consequence of incorporating the contracted Ricci-collineation along the fluid 
flow vector in the Friedmann model (Abdussattar and Vishwakarma 1996; 
Vishwakarma 1999, 2000a). 

The case (2) has been proposed from the similar dimensional arguments 
as made by Chen and Wu (Carvalho et al 1992, Waga 1993) and also from 
other arguments (Lima and Carvalho 1994) and has also been considered 
by several other authors, e.g., Salim and Waga (1993), Arbab and Abdel- 
Rahman (1994), Wetterich (1995), Arbab (1997). In view of the present 
estimates of A being of the order of Hq, this ansatz seems as a natural 
dynamic law for the decay of A. (The subscript zero denotes the value of the 
quantity at the present epoch). 

The dynamical law mentioned in the case (3) has been obtained on the 
dimensional ground by suggesting that A ~ QH 2 which leads to A ~ p, 
(Vishwakarma 2000b), where Q is the density of the matter in units of the 
critical density, i.e., Q = pj p c = 8nGp/3H 2 (note that in our notation, Q 
has contributions from baryons, neutrinos, radiation, etc. and dark matter 
but not from A). 

The case (4) has also been considered by several authors by imposing 
supplementary conditions which are equivalent to assuming a power law for 
the scale factor (Endo and Fukui 1977, Lau 1985, Bertolami 1986, Berman 
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and Som 1990, Berman 1991, Beesham 1994, Lopez and Nanopoulos 1996). 
2. General Features of the models 

In order to find the dynamics of the models resulting from the four cases of 
A mentioned above and then to compare our results with the earlier works, 
we consider the homogeneous and isotropic models in Einstein's theory. We 
assume that the matter source in the universe is described by the perfect 
fluid energy momentum tensor 

T ij = (p + p)u i u j +pg ij , (1) 

where p is the density of baryons, plus neutrinos, plus radiation, etc., plus 
dark matter. The Einstein field equations are, then, given by 

R ij - -Rg ij = -8ttG 
2 y 

For simplicity we consider units with c = 1. In view of the Bianchi identities, 
this equation suggests that A would become constant for T lj = 0. Also, A 
must be a constant if T y is conserved separately (we consider G as a true 
constant). This implies that existence of some non- vanishing matter, how- 
ever small, is necessary to kick A to start varying. In this connection, we 
recall that a decaying A has led to the possibility of non-singular models in 
which initially all the energy is locked up in potential form in the curva- 
ture of spacetime and as the universe evolves, curvature unfolds and creates 
matter (Ozer and Taha 1987, Abdussattar and Vishwakarma 1996). In this 
view, this primordial tiny matter (required by a decaying A) in the other- 
wise empty baby universe might serve the purpose of primordial seeds for 
structure formation. 

We also observe from equation (2) that the conserved quantity is now 
[T y ' — {A(t)/8irG}g t: '] = T t Q tal and not T y and A separately as would have 
been the case, had A remained a true constant. Thus the incorporation of 
a variable A in the Einstein field equations is equivalent to postulating an 
additional source —(A/8irG)g' lj = T % J as the energy momentum tensor of 
vacuum with the energy density of vacuum p v = —p v = A/8nG. 

For the Robertson- Walker metric 

ds > = _ df 2 + 5-2/^ + r 2( d 02 + gin 2 0^2)1 (3) 

1 — kr 2 



rpij 



8nG 



9' 



(2) 
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and the usual barotropic equation of state of the matter source 



p = wp, < w = constant < 1, (4) 

equation (2) yields the following two independent equations: the Friedmann 
equation 

S 2 k 8ttG A 

P+o (5) 



S 2 S 2 3 r 3 
and the Raychaudhuri equation 

S 4ttG A 
-S = — {1 + 3w)p -3- (6) 
It is common in cosmology to consider the constant w — 1/3 for the radiation- 
like case and w = for the dust-like case. Recently —1 < w < has been 
suggested by postulating an 'x-fluid' (Turner and White 1997, Zlatev et al 
1999) to make up the long-sought dark matter. However, these new species 
of matter with negative pressure are equivalent to a variable cosmological 
term and will not be consider here. The divergence of equation (2) leads to 

S 

Ptot + 3(ptot + Ptot) g = 0, (7) 
where p to t = P + Pv and p tot = p — p v . Equation (7) can also be written as 

o-3(l+w) „ 

p = CS~ 3(1+W) — / AS 3{1+w) dt, C = constant. (8) 

8nG J 

This indicates that a decaying A (A < 0) makes a positive contribution to p. 
Equations (5) and (6) can be used to obtain 

2 1 + (1 + 3w) + A) - (1 + w)A = 0, (9) 

which gives the dynamics of the scale factor for a particular A. 

We shall now describe, briefly, the models resulting from the different 
dynamical laws for the decay of A mentioned in the previous section. 
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Case (1): 

We consider 



A = ^, (10) 

where the constant rt\ is a new cosmological parameter replacing the original 
A and is to be determined from the observations. Equations (8), (9) and (10) 
give the dynamics of the model as 

Tli 9~ 2 

p = Cl ^) + -^_ ; , Cl= constant (11) 



and 



£ 2 = %ixGC x s _ (1+3w) ni(3 + «■ 



Case (2): 

We now consider 

A = n 2 # 2 , (13) 

where, as in the previous case, the parameter ti 2 is to be determined from the 
observations. The use of (13) in equations (5) and (8), leads to the following 
equations 



p = c 2 5-(«a-3)(i+t") — n2kS 2 c constant (14) 

1 47rG[(3-n 2 )(l + «;)-2]' V ; 

and 

• 2 87rCC 2 (n2 „ 3)(1+w)+2 _ (l + 3w)fc 

~ (3-n 2 ) D [( 3 -n 2 )(l + «;)-2]' 1 ° J 

which determine the dynamics of the model. 



Case (3): 

We now consider the case 
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A = n 3 p, (16) 

where the parameter n 3 , describing the model, is to be determined from the 
observations as in the earlier cases (1) and (2). 

The dynamics of the models can be obtained from equations (5), (8) and 
(16) as 

p = C 3 S- 3{1+w)/{1+n) , C 3 = constant > (17) 

and 

S 2 = + n ) S^- l ~^ 1+n ) - k, (18) 

3 

where n = ris/SiiG. 



Case (4): 

Finally we consider the case 

A=£, (19) 

where the parameter describes this family of models. By using (19), 
equation (9) can be written as 

x _ 3(1 + wfn A 

At" ' [ Uj 

where the new variable x is related to S by means of the transformation 

S 2 x 



S Z(\ + w)x (21) 

implying that S ~ x 2/3(i+w)_ fj ere we have restricted ourselves to the flat 
model only as has been done by the other authors. By assuming t = e y , 
equation (20) can be further transformed to 

d^x _ dx _ 3(l + w) 2 n 4 ^ _ 

dy 2 dy 4 X " U ' [ZZ) 

which gives a physically viable solution as 
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1 + Jl + 3(1 + w) 2 n 4 

S~t n \ withn 5 = V f- , (23) 

3(1 + w) 

provided n 4 > — 1/3(1 + w) 2 . Other solutions (for n 4 < — 1/3(1 + w) 2 ) are not 
interesting on the ground of age considerations (Overduin and Cooperstock 
1998). Solution (23) is consistent with the different power-law solutions ob- 
tained by various authors by assuming different supplementary conditions, 
as we have mentioned in the previous section. However, in view of H ~ t L 
as is implied by equation (23), this case reduces to the case (2) of A ~ H 2 
and will not be considered separately further. 

While comparing the models with observations, we shall be mostly inter- 
ested in the effects which occurred at redshift z < 5. We, therefore, neglect 
radiation and consider w = in the following. Equations (5) and (6), then, 
can be recast in the following forms to give the relative contributions of 
different cosmological parameters at the present epoch: 

i + n k0 = n + n A0 , (24) 

2[q + n A0 ] = n , (25) 

where Q A = A/3H 2 and = k/S 2 H 2 are, respectively, the energy density 
of vacuum and that associated with curvature in units of the critical density 
and q is the deceleration parameter. 

It would be worth while to find out if there exist some characteristics of the 
models, other than the different dependences of A, which clearly distinguish 
the evolution of the models discussed above. The dynamics of the models, 
described by equations (12), (15) and (18), depend on the parameters rii and 
Cj. In order to make the comparison possible, we write them in terms of H , 
Q an d Q\o as follows: 

ni = 3H 2 S 2 n A0 , C 1 = ^[n -2n A0 ]S 3 0] (26) 

87rG 



n 2 — 3f2 AO , C 2 — — 

87TU 



1 — 3Q A q 



5 3(i-n Ao) . (27) 



n 3 = SvrG^, C 3 = M^f*/^^ (2g) 
i 2q 87rG 
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where So = H 1 yk/Qko with il^o given by equation (24). Equations (12), 
(15) and (18) then reduce, respectively, to the following. 

A ~ S- 2 : 



A~ H 2 : 



A ~ p 



62 tj2 q2 



(too - 2i2 AO )^ + (1 - n + 2O A0 ) 



q2 tj2 q2 

b = N b 



l-3fi A0 \s) l-3fi A0 



q2 tt2 q2 

b = H b 



(n + n A0 ) ( 



g s (Qo-2n A o)/(^0+^Ao) 



+ 1 — Q — Q A0 



(29) 



(30) 



(31) 



It is now clear from equations (29)-(31) that the dynamics of the models 
depend only on the values of Qo and ^ao- Also the relative behaviour of S 
from these equations is quite arbitrary and changes with different pairs of Qq 
and f^Ao- K thus seems that in general there is no feature which distinguishes 
the models. We also note that for f2 = 2fi A (b the models become identical 
reducing to S ~ t and p ~ A ~ t~ 2 . Also, for the trivial case A = 0, the 
models become identical reducing, obviously, to the standard FRW models 
with zero A. 

We shall now derive the m-z and Q-z relations in the models in order to 
compare them with observations, which will be done in the following sections. 



3. Study of magnitude-redshift relation in the models 

Suppose that the observer at r = and t — t receives the light emitted at 
t — ti from a source of absolute luminosity L located at a radial distance 
r\. The cosmo logical redshift z of the source is related with t\ and t by 
1 + z — S(t ) I S(ti). If the (apparent) luminosity of the source measured by 
the observer is /, the luminosity distance of the source, defined by 

is 4^p (32) 

is then given by 
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d L = (1 + z)S n. (33) 

For historical reasons, the absolute and apparent luminosities L and I are 
defined, respectively, in terms of the absolute and apparent magnitudes M 
and m as L = 1(T 2M / 5 x 3.02 x 10 35 erg s" 1 and I = lO" 2 " 1 / 5 x 2.52 x 10" 5 erg 
cm -2 s _1 (Weinberg 1972). When written in terms of M and m, equation 
(32) yields 



m(z] M, too, Q A0 ) = M + 51og 10 



V L (z; fi , ^ao) 



(34) 



IMpc 

where M = M — 51og 10 i7 + 25 and ^o? ^ao) = #o ^l(z; ^ao? #o) is 

the dimensionless luminosity distance. By using equation (3), the coordinate 
distance r 1 , appearing in equation (33), yields 

^)=fi n M (35) 

JS /(l+z) SS 

with 

tp(ri) = sin _1 ri, k — 1 
= n, A; = 

= sinh~ 1 r 1 , k = —1. (36) 

The luminosity distances for the different models, described for the cases 
(1), (2) and (3), will be given by calculating the right hand side of equation 
(35) for these models. By using equations (29)-(31), equation (35) can be 
written, for the three models, as 

^i) = ^V J .> ( 37 ) 

where Ij (j — 1, 2 and 3 correspond to, respectively, the three models de- 
scribed for the cases (1), (2) and (3)) is given by 

h = f [(2ft A0 - Q + 1)(1 + z'f - (2fi A0 - n )(l + z'fY^dz', (38) 
Jo 
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'(Q - 2VL 



AO ) 



3Q 



AO) 



(1 + z'f 1 -^ + 



(1 - 3fi A0 ) 



AO J 



1-1/2 



(1 



J\2 



dz! 
(39) 



and 



h= [(fi + fi AO )(l + /) 3 ^0+^«) + (l_ fio _ fiAO )(l + ^)2 r l/2^ (4Q) 
JO 

(For convenience, we also write the similar integral I c for the constant A- 
Friedmann models obtained as 

I c = f [(1 + Q z')(l + z'f - z'(2 + z'^-^dz', 
Jo 

which will be used to compare our results with those for the constant A- 
models.) 

Equations (33), (36) and (37) can also be combined into a single compact 
equation to give 

V L (z; flo, fi A0 ) = £ (VKIj) , (41) 

where 

£(x) = sin(x) with /C = Qko when Q k0 > 0, 

£(x) = sinh(:r) with /C = — Qko when fl^o < and 

£(x) = x with K, — 1 when f2fc = 0. 

Thus for given Ai, and £7ao, equations (34) and (41) give the pre- 
dicted value of m(z) at a given z. Using the observed values of the effective 
magnitude rrf^ (corrected for the width-luminosity relation) and the same 
standard errors a z ; L and a e s of the ith supernova with redshift Zi as used 

by Perlmutter et al, we compute x 2 according to 

, " [mf-m{ Zl )f 



The best fit parameters are obtained by minimizing this equation. 

We consider the data set on m eff and z of 54 supernovae as used by 
Perlmutter et al in their primary fit C. This set comprises 38 high-redshift 
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supernovae from the Supernova Cosmology Project (excluding 2 outliers and 
2 likely reddened ones from the full sample of 42 supernovae) together with 16 
low-redshift supernovae from the Calan-Tololo sample (excluding 2 outliers 
from the full sample of 18 supernovae). We fit the low- and high-redshift 
supernovae simultaneously to equation (34) to estimate the parameters At, 
Q and f^Ao- This fitting procedure is followed by Perlmutter et al in their 
primary fit C. We have listed the best-fit values in Table 1. The global best- 
fit parameters have been calculated by giving free rein to Qo, Q\o an d At, 
i.e., at 54 — 3 = 51 degrees of freedom (dof) whereas the best-fit parameters 
for flat models are calculated at 52 degrees of freedom with the constraint 
Qo + ^ao — 1- For comparison, we have shown the best-fit parameters 
for the constant A-models. The x 2 values obtained at different degrees of 
freedom have also been converted into confidence levels (cl) [the goodness of 
fit probability Q = 1 — cl%] . In Figure 1, we have shown the fit of the data 
to the best-fitting flat models and compared it with the fit of the data to the 
Einstein-de Sitter model. 

It may be mentioned that Perlmutter et al have also fitted the data for 
only 3 parameters Ai, flo and Aao an d not for the 4 parameters Ai, flo, ^ao 
and a (the slope of the width-luminosity relation) as mentioned in their paper 
(from a personal discussion with Professor R. S. Ellis, one of the authors of 
the paper). A self-consistent 4-parameter fit has been done by Efstathiou et 
al (1999). 

Table 1. Best-fit values of different parameters for three variable A-models 
estimated from the Type la supernovae data. For comparison, the best- 
fit parameters for the constant A-Friedmann model (Perlmutter et al') and 
Einstein-de Sitter model have also been shown. 



Models Best-Fit Parameters (Flat models) Global Best-Fit Parameters 
fl M x 2 dof cl(%) f] n A0 M x 2 dof cl(%) 



A~ S- 2 


0.49 


23.97 


58.97 


52 


76 


1.86 


1.52 


23.95 


57.43 


51 


75 


A~ H 2 


0.4 


23.96 


58.33 


52 


75 


0.98 


1.53 


23.91 


56.86 


51 


73 


A ~ p 


0.4 


23.96 


58.33 


52 


75 


1.62 


1.59 


23.93 


56.91 


51 


74 


A = const 


0.28 


23.94 


57.68 


52 


73 


0.79 


1.41 


23.91 


56.82 


51 


73 


A = 0, 

n = i 




24.21 


92.89 


53 


99.9 
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Figure 1. Hubble diagrams for 38 high-redshift and 16 low-redshift su- 
pernovae are shown. The solid curve represents the best-fitting flat model 
Qo = 0.49 (for A ~ S~ 2 ) and the dotted curve represents the best-fitting flat 
model f2 = 0.4 (for A ~ H 2 ). The curves overlap on each other because of 
the same goodness of fit. For comparison, the canonical Einstein-de Sitter 
model (Vl k = 0, ^a = 0) has also been plotted (dashed curve). 

The fits of the data to all the three variable A-models are almost as good 
as to the constant A-Friedmann model (Perlmutter et al') with x 2 /dof ~ 1 
(a very good fit indeed) and require non-zero, positive values of A. It is clear 
from Table 1 and Figure 1 that the data do not fit the Einstein-de Sitter 
model (f2 = 1,A = 0). We note that the global best-fit values for f2 and 
f^Ao are always larger than those for the corresponding flat models. Also 
the estimates of Qq for the variable A-models are always higher than that 
for the constant A-Friedmann model (Perlmutter et al') and hence require 
dark matter, though not as much as in the non-flat models. It is also clear 
from the Table that for all the three variable A-models, f2 /2 is less than the 
corresponding ^ao (as is the case for the constant A-Perlmutter et al' model) 
leading to q < 0. 
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One may notice that the estimates of the parameters for the flat models 
for the cases A ~ H 2 and A ~ p are exactly same. This simply reveals the 
fact that for k = 0, the two models become identical. This can be checked 
from equations (17) and (18) which, for k = 0, give S ~ f2(i+n)/3(i+w) \ ea dmg 
to A ~ p ~ t~ 2 ~ H 2 . This may be considered as a cross check for our fitting 
procedure. 

Elsewhere (Vishwakarma 2000a) we have used this data set with a dif- 
ferent fitting procedure for the model A ~ S~ 2 . Instead of fitting the high 
and low-redshift data simultaneously, we used low-redshift measurements to 
determine Ai and by using this value in equation (34) we performed a two- 
parameter fit of the high-redshift supernovae to estimate Qq and f^Ao- The 
best-fit values so obtained are in good agreement with the estimates shown 
in Table 1 for that model. 

We would like to mention that the supernovae data have also been used 
in many studies to place constraints on models with 'x-fluid' (p x = wp x 
with — 1 < w < 0) or 'quintessence' (a slowly varying scalar field) (Waga 
and Miceli 1998; Efstathiou 1999; Perlmutter, Turner and White 1999; Po- 
dariu and Ratra 1999; and the references therein). Though 'x-fluid' and 
'quintessence' are formally equivalent to a variable A, however, these models 
are essentially different from our kinamatical 'usual' A-models considered in 
this paper in the sense that in our models matter is created as a result of 
the decaying A, as has been mentioned in section 2, whereas in the exotic 
fluid description, the exotic component is conserved. In the present approach 
(or in the case of 'x-fluid' with w = —1), if the A-source [{A(t)/8irG}gij] is 
conserved separately, then A would become constant, as has been mentioned 
earlier. 

4. Study of angular size-redshift relation in the models 

The compact radio sources, of angular size of the order of a few milliarcsec- 
onds, which are observed with very long baseline interferometry (VLBI), are 
deeply embedded in the galactic nuclei and are expected to be free from the 
evolutionary effects. These sources are attractive for the test of cosmological 
models because they are usually identified with quasars which, in general, 
have high redshifts, so the differences between models may be more easily 
distinguished than for the extended double-lobed sources. 

Kellermann (1993) used a sample of 79 such sources and showed that a 
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credible Q-z relation can emerge. His main motivation was, however, limited 
to showing that the resulting Q-z diagram was cosmologically credible and 
compatible with the Einstein-de Sitter canonical model with A = rather 
than to be specific about the acceptable ranges for the various cosmological 
parameters. Extending his work, by introducing A into the picture, Jackson 
and Dodgson (1996) showed that the Kellermann's data could also fit very 
well with a low density, highly decelerating model having a large negative 
cosmological constant. Later on, a more extensive exercise was carried out 
by Jackson and Dodgson (1997) by considering a bigger sample of 256 ultra- 
compact sources selected from the compilation from Gurvits (1994). Then, 
they concluded that the canonical cold dark matter model with A = was 
ruled out by the observed relationship at 98.5 percent level of confidence and 
low density FLRW models with either sign of A were favoured. Elsewhere we 
(Vishwakarma 2000b) used the data set from Jackson and Dodgson (1997) in 
the model resulting from the case (3) and found the model in good agreement 
with the data. 

Recently this data set has been updated and extended by Gurvits et al 
(1999) as a sample of 330 sources (though at the expense of homogeneity) 
distributed over a wide range of redshifts 0.011 < z < 4.72. In order to 
minimize any possible dependence of linear size on luminosity and that of 
angular size on spectral index, they discarded lower values of luminosities and 
extreme values of spectral indices and selected only 145 sources with total 
radio luminosity Lh 2 > 10 26 W/Hz and —0.38 < spectral index < 0.18, which 
are hoped to be free from the evolutionary effects and hence conceivably 
comprise a set of standard objects. This sub-sample was distributed into 
12 bins, each bin containing 12-13 sources. They studied the Q-z relation in 
Friedmann models with A = for this sample and found the best-fit model as 
go = 0.21 with no-evolution assumption. However, as we shall see, the models 
with a non-zero A (constant) are also in good agreement (even better) with 
the data. We shall also use this data set to study the Q-z relation in the 
variable-A models considered in this paper. For this purpose, we now derive 
the Q-z relation in the models by assuming that there are no evolutionary 
effects due to a linear size-luminosity or linear size-redshift dependence. 

The angular diameter distance d& of a light source of linear size d and 
angular size 6, defined by 
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~ 0' ^ 

is related to the luminosity distance of the source by 

By comparing this with equation (33) and using (36), (37) and (43), we 
obtain the Q-z relation in the models as 



e(z; dh, Q Q , n A0 ) = (lOO km s^Mpc" 1 ) dhVK, (1 + z) / £(VlC Ij), (45) 

where h is the present value of the Hubble parameter in units of 100 km s _1 
Mpc -1 . 

It is thus clear that for given dh, Cl an d ^ao, the predicted value Q(z), 
at a given z, is completely determined. Using median values of the grouped 
data into 12 bins, as used by Gurvits et al, we calculate the theoretical Q(z) 
for a wide range of parameters dh, Qo and Qao and compute x 2 according to 

2 " r e, - &( Zl ) 

x hi * 

which measures the agreement between the theoretical Q(z) and the observed 
value @j with the errors <jj of the ith bin in the sample. The maximum 
likelihood estimates of the parameters are obtained by minimizing x 2 as 
usual. We have listed our results in Table 2 below. 



(46) 
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Table 2. Best-fit parameters for the Friedmann model (with constant A) and 
three variable A-models estimated from Gurvits et al' data. For comparison, 
the best-fit parameters for Gurvits et al' model have also been shown. (The 
values of dh are shown in unit of 1 pc.) 



Models Best-Fit Parameters (Flat models) Global Best-Fit Parameters 

tt dh x 2 dof cl(%) tt Vl A0 dh x 2 dof cl(%) 



A~S- 2 


0.68 


18.25 


4.71 


10 


9 


0.97 


0.61 


17.99 


4.66 


9 


14 


A ~ H 2 


0.67 


18.53 


4.72 


10 


9 


0.29 


1.03 


23.71 


4.28 


9 


11 


A ~ p 


0.67 


18.53 


4.72 


10 


9 


0.53 


0.82 


22.41 


4.27 


9 


11 


A = const 


0.2 


22.29 


4.52 


10 


8 


0.08 


1.16 


27.73 


4.17 


9 


10 


Gurvits's 
model 












0.43 





17.72 


4.74 


10 


9 



We thus find that the Gurvits et al' model (Friedmann model with a van- 
ishing A) is not the best-fitting model. The best-fitting model with constant 
A is a low-density model dominated by positive A. It is clear from Table 
2 that all the models, variable- or constant-A, show an excellent fit to the 
data with x 2 /dof ps 0.5 and require non-zero, positive values of A. Here 
again the best-fit values of f2 for variable A-models are larger than the best- 
fit f2 for the constant A-model. Unlike the supernovae data, this data set 
requires both, the accelerating and decelerating models: the global best-fit 
parameters (for the variable as well as the constant A-models) represent an 
accelerating expansion whereas the best-fit parameters for the flat models 
(only the variable A-ones) represent a decelerating expansion. 

Once again we encounter with the same values of the parameters for the 
models A ~ H 2 and A ~ p for the flat cases, which serves a cross check for 
the fitting procedure. In Figure 2, we have shown the fit of the data to the 
flat models. 
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Figure 2. Gurvits et al' data points are compared with the three best-fitting 
flat models obtained for the cases: A ~ S~ 2 with Q = 0.68 (solid curve), 
A ~ H 2 with Q = 0.67 (dotted curve) and A = constant with Q = 0.2 
(dashed curve). 

Following Gurvits et al, it would be worth while to study the consequences 
if the source linear size depends on the source luminosity and redshift (i.e., 
the sources are not the 'true' standard rods). For this purpose, we consider 
the following phenomenological expression for as has been considered by 
Gurvits et al: 

= d/d A ~ £h 1/(1 + zyV A \ (47) 

where Va = -f^o^A is the dimensionless angular diameter distance and £h is 
the linear size scaling factor which reduces to dh for the no-evolution case. 
The parameters (3 and 7, representing the dependence of the linear size on 
the source luminosity and redshift respectively, are to be estimated from the 
data. We fit the Gurvits et al' sub-sample to equation (47) for five parameters 
Q , ^h, j3 and 7. The best-fit values for the three models are obtained 
as follows. 
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A ~ S- 2 : 

tt = 0.89, tt A0 = 0.58, eh = 23.69, (3 = 0.37, 7 = -0.45, x 2 = 3.21; 



A~H 2 : 

n = 0.39, n A0 = 0.95, eh = 28.16, p = 0.37, 7 = -0.44, x 2 = 2.72; 

A ~ p : 

fi = 0.51, fi A0 = 0.83, eh = 17.9, /3 = 0.38, 7 = -0.42, x 2 = 2.47. 

We find that the estimates of f2 and £Ia thus obtained are in agreement 
with those shown in Table 2. The estimates of (3 and 7 are in qualitative 
agreement with (5 = 0.37 and 7 = —0.58 obtained by Gurvits et al. The 
values of Qq and Q^o for the case A ~ p are also in qualitative agreement with 
the estimates of these parameters obtained (by assuming no-evolution) from 
an independent sample of sources (Vishwakarma 2000b). It thus seems that 
there is no need to consider evolutionary effects due to linear size-luminosity 
or linear size-redshift dependence. However, this result, which is based on 
the above-made comparisons, is rather superficial than indicative because 
there are large uncertainties in the above-mentioned estimates which are due 
to the large number of regression parameters. 

5. Conclusion 

We find that all the variable A-models, fit the data sets equally well and a 
discrimination between the different dynamical laws or between the different 
geometries is not possible with the present status of the data sets. While the 
fit of the models (constant, as well as, variable A) to the supernovae data 
is very good with x 2 /dof ~ 1, their fit to the radio sources data is excellent 
with x 2 /dof even less than I. The estimates of the density parameter Q 
(from both data sets) for the variable A-models are found higher than those 
for the constant A-Friedmann model. The supernovae data set requires an 
accelerating expansion for the variable A-models as it does for the constant 
A-model. However, the radio sources data set requires accelerating as well as 
decelerating expansion for the variable A-models, though only an accelerating 
expansion for the constant A-model. 

Many of the fl values shown in Tables 1 and 2 (especially the global best- 
fit estimates from the supernovae data) are found rather high. In view of the 
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small observed value of Qq, this implies that not all the best-fit solutions 
found here are realistic. For 0.2 <~ Qq <~ 0.5, which seems to be the 
general consensus at present, the only realistic solutions are (a) estimated 
from the supernovae data: the best-fit solutions for flat models (including 
the constant A-case); (b) estimated from the radio sources data: the global 
best-fit solutions for the models A ~ H 2 and A ~ p, the best-fit solution for 
the flat model with A = constant and the Gurvits et al' model. 

As in the case of the recent analyses of the CMB results from BOOMERANG 
and MAXIMA (Jaffe et al, 2000), our global best fit values from both data 
sets imply a spherical universe (Qq + ^ao > 1). 

One may notice from Tables 1 and 2 that the best-fit values of f^o and f^Ao 
from the two data sets are a bit too different for the corresponding models, 
especially in the non-flat cases. The reason for this may be the incomplete 
understanding of a number of astrophysical effects and processes (evolution, 
inter galactic dust, etc.). It should also be noted that the Gurvits et al' com- 
pilation is based on very inhomogeneous data obtained by many different 
observers using different instruments and imaging techniques. New VLBI 
observations now in progress will improve the accuracy of the observed O — z 
relations as it will provide a uniform data set for analysis. A reasonable hope 
is that, as the data accumulates and its understanding improves, such tests 
coupled with plausible evolutionary hypothesis will lead to similar results. 
These results should be consistent with the CMB measurements and also 
with the spatial scale of the peak in the power spectrum of galaxy pertur- 
bations which provides a standard ruler in co-moving space for measuring 
cosmological parameters (Broadhurst and Jaffe 2000, Roukema and Mamon 
2000). 
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